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Abstract. The traditional perturbative method is ap- 
plied to the case of gravitational lensing of planetary sys- 
tems. A complete and detailed description of the structure 
of caustics for a system with an arbitrary number of plan- 
ets can be obtained. I have also found precise analytical 
expressions for microlensing light curves perturbed by the 
presence of planets. 
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1. Introduction 

In the last years, the prediction and the observation of 
many microlensing events are gathering ever more inter- 
est in gravitational lensing. The typical light amplifica- 
tion curve for these events, found by Paczynski in 1986, 
has been observed by several astronomical collaborations 
in observation campaigns toward the bulge of our galaxy 
(Udalski et al., 1993a; Alard et al., 1997), the Large Mag- 
ellanic Cloud and the Small Magellanic Cloud (Alcock 
et al., 1993; Aubourg et al., 1993), the spiral arms and 
Andromeda galaxy (Tomaney & Crotts, 1996; Ansari et 
al., 1997; Melchior et al., 1999). Recently, observations 
toward globular clusters have even been suggested (Jet- 
zer, Strassle & Wandeler, 1998). Beyond proving the cor- 
rectness of Paczynski's predictions, the observation of mi- 
crolensing events provides a very cunning instrument for 
the investigation of the halo of our (and/or some other) 
galaxy. 

Together with events strictly following Paczynski's 
curve, some events showing deviations from the standard 
behaviour have been detected. Each of these deviations 
has found some interpretation (Finite source cut-off (Witt 
& Mao, 1994; Alcock et al., 1998), blending (Sutherland, 
1998), parallax effect (Gould, 1992; Alcock et al., 1995), 
binary lens (Schneider & Weiss, 1986; Mao & Paczynski, 
1991; Udalski et al., 1993b). Indeed, the most intriguing 
of these deviations is the one induced by a binary (or mul- 
tiple) lens. 
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The study of light amplification curves produced by 
multiple lenses has not yet been performed analytically 
because of the difficulties in the inversion of the lens ap- 
plication. Anyway, these curves can be obtained numeri- 
cally by using some inversion algorithm (inverse ray shoot- 
ing) (Wambsganns, 1997). From the analytical point of 
view, only the caustics of a general binary lens have been 
studied in some detail (Witt & Fetters, 1993). The lack of 
analytical results utilizable in microlensing constitutes an 
irksome obstacle in the complete interpretation of multiple 
microlensing events. 

A particularly interesting case of multiple 
Schwarzschild lenses is formed when one of the masses 
is much biggest than others (Mao & Paczynski, 1991; 
Gould & Loeb, 1992). This is the situation of a typical 
planetary system where a central star is surrounded by 
its planets bearing masses thousand or million times 
smaller. The perturbations on Paczynski's curve induced 
by the presence of a (even Earth - like) planet are in 
principle detectable by collaboration teams exploiting 
world wide telescopic networks (Peale, 1997; Sackett, 
1997). Then microlensing could become a new efficient 
method for the detection of small planets in extra - solar 
systems. This justifies the major interest in this field that 
is growing in the last months. Preliminary calculations 
on the probability of detection of planets have been 
made (Gould & Loeb, 1992; Bolatto & Falco, 1994) and 
great efforts are lavished on the problem of extraction of 
planetary parameters by approximate models (Gaudi & 
Gould, 1996). It is easy to imagine how the availability 
of an analytical expression for light curves could help the 
researches in this field. 

The aim of this work is to describe planetary ef- 
fects perturbatively, exploiting the very little ratios of the 
masses of the planets with respect to the star mass (Gould 
& Loeb in 1992 first pioneered this kind of approach). In 
section 2 the lens equation and other usual objects are 
specified for the case of planetary systems. In section 3, by 
means of perturbative theory, I derive the complete struc- 
ture of critical curves and caustics of a general (not only 
binary) planetary system; position of planetary caustics 
and cusps are also found. In section 4 the problem of the 
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Fig. 1. Tipical scheme of planetary lensing. 



inversion of lens application is faced and resolved; conse- 
quently, analytical microlcnsing light curves for planetary 
events thus obtained are shown. 



2. Planetary lensing 

A planetary system is nothing but a discrete set of 
Schwarzschild lenses in a small portion of the space. Fig. ^ 
shows the typical situation of planetary lensing. The lens 
plane is defined as a plane orthogonal to the line of sight 
situated on the barycentre of the lens mass distribution. 
According to the ordinary theory of lensing (Schneider, 
Ehlers & Falco, 1992), if the scale of this distribution is 
much smaller than the distances separating the lens from 
the source and the observer, then one can deal with the 
density projected on the lens plane instead of consider- 
ing the original volume density. This hypothesis is quite 
verified in real observations where the typical distances 
are at least of the order of kpc. This consideration can 
be important in planetary systems where very far plan- 
ets can eventually have projections enough near to the 
central star to give perturbations comparable to those of 
planets placed in more favourable positions. This means 
that multiple events could be less out of common than one 
could think (Wambsganns, 1997; Gaudi, Naber & Sackett, 
1998). So, it is desirable to preserve the whole planetary 
system as much as possible before abandoning it for the 
simplest case of the single planet around the big star. We 
shall see that the perturbative theory has the consider- 
able advantage of being rather insensitive to the number 
of planets as regards the difiiculty of the problem. 
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Let's define the length — ^, — 2 rr 

{xi]X2) shall denote the coordinates in the lens plane nor- 
malized to i?^, while y = (?/i; 2/2) shall be the coordinates 
in the source plane normalized to i?^-^^. mi will be the 
mass of the central star and 7712, m„ will be the masses 
of the planets. All of these masses are meant to be mea- 
sured in solar masses. The star will always be placed at 
the origin, while the projection on the lens plane of the 
ith planet will be denoted by Xi = {xii;xi2). With these 
notations, the lens equation reads: 



mix 



E 



rui (x - Xi) 



(1) 



In (|]) , the deviation of light rays due to the star has 
been explicitly separated by those of the planets. Given 
a source position y, the corresponding x solving the lens 
equation are called images. 

Many interesting properties of this vectorial applica- 
tion can be studied through its Jacobian matrix. In par- 
ticular, the determinant of this matrix contains nearly all 
the information about the properties of the images. Let's 
write the Jacobian determinant for the case of planetary 
lensing using ffl): 



det J = 1 
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where Ai = (A^i; Ai2) = x — x^. Given an image I at 
position X/, the sign of det J (x/) is called the parity of I. 
It can be proved that the amplification of the image I is 
given by 



1 



|det J(x/) 



(3) 



We see from this equation that when det J is null, the 
amplification diverges. This is rigorously true for point 
sources (in ray optics), while for finite (real) sources the in- 
tegration over the source's surface makes the amplification 
finite (Witt & Mao, 1994). The points where det J van- 
ishes are called critical points, the corresponding points 
in the source plane through (|^) constitute the caustics. 
So, a point source crossing a caustic will produce images 
with infinite amplification. Real sources crossing caustics 
by some of their parts will be highly (but not infinitely) 
amplified. 

The structure of critical curves and caustics provides 
a substantial description of the general behaviour of the 
lens (number and roughly location of the images can be 
established). In microlensing, it is possible to give a qual- 
itative description of light curves just checking whether 
the track of the source threads some caustic or not. 
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In fig. Q the (numerically obtained) critical curves and 
the caustics of a star with a single planet are shown. There 
are three possible cases in such a situation. I recall that 
for a single point source the critical curve is a ring with 
radius given by its Einstein radius \/m^ while the caustic is 
a point in the origin of the source plane. When the planet 
is far beyond the star's Einstein ring, there is only a small 
perturbation of the two rings which lends finite extension 
to the originally point - like caustics. This is much more 
evident in the planetary caustic which is also displaced 
towards the star. When the planet is in the proximity of 
the star's Einstein ring, the two critical curves merge and 
so do the caustics. In the last situation where the planet is 
internal to the star's Einstein ring, the star's critical curve 
returns to be very near to a ring while the planet's critical 
curve turns into two ovals to which a couple of triangular 
caustics correspond behind the star. 





3. Caustics and perturbative analysis 

In the solar system, the mass ratios between planets and 
the sun are always less than one thousandth. Jupiter 
is 9.4 X WT'^Mq. Other planets are even less: Earth is 
3 X 1O^^M0. With these numbers, it is natural to ex- 
pect that the presence of planets should cause only little 
perturbations to the single lens case. Upon this considera- 
tion the perturbative hypothesis is based. In this and the 
following section the ratios between planetary and stel- 
lar masses will play the role of perturbative expansions 
parameters. We shall see that in most cases a first order 
expansion is sufficient to get very reliable results. 

Let's turn to the caustics of planetary systems. First 
I shall examine the modifications induced in the star's 
Einstein ring and consequently the central caustic. Then 
I shall deal with planetary caustics. 
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Fig. 2. There are three possible situations in planetary 
lensing as regards critical curves, shown in these plots. 
The smooth round curves are critical curves, while the 



3.1. Central caustic 

Of course, the starting point for the study of critical curves 
is the equation det J = 0, which can easily be written in 
polar coordinates: 



mi (cos^ 7? 



+ 4 



mi sini^cosi^ 



= 1 (4) 



Here Ai — (r cos — xn ; r sin ■& — Xi2). Expanding this 
equation to the first order in m^, we get: 
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1 (5) 
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Fig. 3. Central caustic for Jovian 
planets. The caustics on the left are 
those found perturbatively, while 
those on the right are obtained nu- 
merically. Little differences occur 
when one of the planets is close to 
the star's Einstein ring. 



The zeroth order solution is simply r — y/mi, i.e. the 
Einstein ring. Let's write the first order solution as: 

r = ^{l + e{d)) (6) 

with e <C 1. Substituting in (|^) and expanding to the first 
order in e. The zeroth order solution cancels and e is found 
solving the remaining first degree equation: 



e (^) = 1 (cos^ d - sin^ d) ^ 



" m AO AO 

2sin^cos,?^^-^^^pi^=i^ (7) 



i=2 



[(A°i)V(A°2) 



where = (y^roTcosi? — 2:^1; y^mi sin?? — 2:^2 ) is to 
the zeroth order. 

By very few steps we have found the perturbation of 
the Einstein ring in a very simple way. The parametric 



equation of the central caustic is soon found by applying 
the lens equation (|^) and expanding again to the first order 
in m, : 
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2 [(AO,f + (A0, 



(8b) 



Of course, perturbative results are characterized by 
precise limits of validity. In our case we see that when 
the denominators in (0) vanish the perturbation diverges. 
This is not allowed by our assumption that e must be very 
small with respect to the unity. Those denominators repre- 
sent the distance between the planet and the general point 
of the unperturbed star's Einstein ring. So we expect the 
perturbative theory to fail in those portions of the critical 
curve which are very near to one of the planets at least. 



Valeric Bozza: Perturbative analysis in planetary gravitational lensing 



5 



We can understand this "failure" if we look back at fig. ||b: 
when the planet is close to the star's Einstein ring, there 
is only one critical curve which is somewhat different from 
the ring in the proximity of the planet. For some values 
of "(5, the radial coordinate describing the critical curve 
assumes also more than one value; this situation cannot 
be described by a first order approximation, where, as we 
saw, the perturbation solves a first degree equation. 

The most interesting aspect of eqs. (0) and (||) is that 
they are comprehensive of the action of the whole plan- 
etary system: they are valid for an arbitrary number of 
planets, not only the classically investigated case of the 
single planet. So these formulas enjoy a very high gener- 
ality and can be used in more realistic contexts. We also 
note that the contributions coming from different planets 
superpose without interfering. This is an obvious conse- 
quence of the first order approximation; if I had included 
second order terms, I would have found "interaction" be- 
tween planets. These "interaction" terms are thus not rel- 
evant in a first approximation. 

Now let's compare the perturbative caustics with those 
found by classical numerical algorithms to test the validity 
of the perturbative approach. In fig. || I show the results 
for the case of two Jovian planets placed in several posi- 
tions. When the planets are far enough from the Einstein 
ring (fig. 1^, ||d), the caustic found according to (||) is en- 
tirely identical to the numerical one. Letting one of the 
planets approach the Einstein ring, a small deviation be- 
gins appearing in the region coming from the portion of 
the Einstein ring that is closest to the planet. This devi- 
ation manifests itself in the size of the largest cusp. For 
Jovian planets these discrepancies unveil at distances from 
the star's Einstein ring of the order of a tenth of the Ein- 
stein radius. These first encouraging results become much 
better in the case of Earth - like planets. We expect the 
range of validity of perturbative results to be increased for 
this kind of planets, because of their smaller mass. Fig. || 
shows that for these little planets things go very well down 
to a hundredth of the Einstein radius. 

So the perturbative method is likely to provide reliable 
results at the first order already. Moreover, it is not to be 
forgotten that, in principle, the approximations can be 
improved pushing farther the perturbative expansion. 

3.2. Planetary caustics 

Planetary caustics are usually studied considering the 
planet as a point-lens with an external shear due to the 
star's gravitational field (Schneider, Ehlers & Falco, 1992). 
This kind of lens was introduced by Chang & Refsdal 
(1979; 1984) in a cosmological context. However, this lens 
is valid in planetary systems only to the lowest non-trivial 
order in m^. Therefore a correct study should only retain 
the lowest order terms in the critical curves equation, so 
that Chang & Refsdal's caustics are a suitable approxi- 
mation at the lowest order only and not beyond. 




Fig. 4. Central caustics for Earth - like planets. The per- 
turbative ones are on the left and the numerical on the 
right. 



In order to complete the discussion of caustics in plan- 
etary systems and study their features properly, in this 
subsection I derive planetary caustics from perturbative 
hypothesis paying full attention to the order of each term. 
The situation for planetary caustics is rather different 
from that of the central caustic. There is no zeroth order 
solution to start from, since their very presence is pertur- 
bative. Nevertheless this is not a great problem: in fact 
we shall just search for the lowest order solution of the 
critical curves equation. 

To achieve this, I now rewrite det J = in polar co- 
ordinates choosing the planet TO2 situated in X2 as the 
origin: 



6 



Valerio Bozza: Perturbative analysis in planetary gravitational lensing 



7712 (cos^ •& — sin^ 1?) 



E 



Af^)' 



7712 sin i5 cos d 
+ n 



"iiAiiAi2 



(A?i 
miAiiAi2 



A?2) 



AI2) 



that what counts is actuaUy the projection of the positions 
on the lens plane. So planets could be very far apart but 
have near projections yielding exotic critical curves. 

We see that the critical curves traced by (|ll]) have two 
branches according to the double sign in their expression. 
For planets external to the star's Einstein ring { Xp > 
^mi), the branch coming from the positive sign is real 
while the other coming from the negative sign is imaginary 



When the planet is very far from the star, we know 
that its critical curve tends to an Einstein ring with ra- 
dius r = Y^mJ. So we search for critical curves solving (J^ 
with r ^ y/m2. Let's save the lowest order terms only. In 
this operation, the contributions coming from the other 
planets are ejected out from the equation. It is convenient 
to place the star in the usual position {—Xp; 0). What re- 
mains is: 



— 1 (9) for all values of 1?, being 
For internal planets ( x. 



2mim2 



(cos^ d 



sm' 



7?) = 1 



(10) 



which is biquadratic in r. The solution is: 



TO2- 



cos2z9± Jl - ^sin^ 2^? 



1 - 



(11) 



which verifies our assumption r ^rni- The parametric 
equations of caustics can be found in the usual way sub- 
stituting ( |ll| ) in the lens application and expanding to the 
first non trivial order (^7712): 



"ii / r 
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mi 



r sin i) 



TO2\ 



sint? 



(12a) 
(12b) 



The contributions from the other planets are again of 
higher order. So the structure of planetary caustics is not 
affected by the presence of other planets at the lowest or- 
der in a perturbative expansion. These formulas can thus 
be used in a single planet situation as well as in a rich 
planetary system. 

Observe that r goes to infinity as Xp tends to ^/rr^l, i.e. 
when the planet is next to the star's Einstein radius. So 
the perturbative results will not be valid in this situation. 
The reason is the same discussed for the central caustic. 
The merging of the two caustic is not describable in the 
lowest order perturbative expansion. Moreover, there's an- 
other limit to be taken in account. I have eliminated all 
the terms coming from the other planets because of their 
higher order. But these terms can become dominant when 
their denominators are small. This happens when one of 
these planets is close to the planet we are examining. This 
is not an exotic situation since we must always remember 



negative and 



^ cos 27? 




< Wl - ^sin^27?. 
the denominator is 



sm 



^ 27?. So the two 



branches are both real for: 



sin^ 27? < 



cos 27? < 



(13) 



that is in two small regions around 7? = ^ and 7? = 
They are both imaginary elsewhere. All these results are 
coherent with the behaviour exposed in fig. ||. We have 
one planetary caustic for external planets and two discon- 
nected caustics for internal planets. 

Fig. 1^ shows the comparison with the numerical caus- 
tics. In fig. 1^ the discrepancy with the numerical results 
appears as a loss of symmetry of the numerical caustic 
which is elongated towards the central star. This effect is 
not present in the perturbative Chang & Refsdal's one. For 
internal planets near the star's Einstein ring, the basis of 
the triangular perturbative caustic is parallel to the star - 
planet axis ( fig. |5|:). So Chang & Refsdal's lens works good 
until the field can be taken as uniform. When the spher- 
ical symmetry becomes important, the caustics begin to 
differ from perturbative ones. These effects can be taken 
into account by considering higher order terms in the ex- 
pansion. These terms would provide the right corrections 
to the Chang & Refsdal's approximation. 

Eqs. (^ij) and (^2|) can be employed to find interest- 
ing characteristics of planetary caustics. For example, let's 
find the position of the couple of caustics for internal plan- 
ets. We saw that the critical curves are centered upon 
7? = ^ and 7? = Consider the first of these (the other 
is similar at all). Inserting these values of 7? in (pi]), the 
possible values of r are: 



\ 



7772 



(14) 



The point 



777 i 7772 



771T 



X„ 



(15) 



obtainable by a quadratic mean from the two values, is 
internal to the critical curve and gives an approximation 
for its position. Immediately, using the lens equation and 
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Fig. 5. Planetary caustics for a Jo- 
vian planet. The caustics on the left 
column are perturbative while the 
ones on the right are numerical. 



expanding to the lowest order terms, we find the position 3.3. Cusps 



of the caustics 

yi 



y2 



Up 



nil 

X. 
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miX: 
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(16a) 
(16b) 



The first of these is a well known formula (Griest & 
Safizadeh, 1997). The second completes the information 
given from the first and allows a complete individuation 
of the two caustics. Fig. ^ is a plot of the position of the 
two caustics as a function of the distance of the planet 
from the star. When Xp —^ the caustics approximately 
move on the lines: 



y2 



V "ii 



(17) 



The two planetary caustics delimitate a region of high 
de - amplification which can appear in microlensing light 
curves as negative peaks. The positions of the caustics can 
give a measure of the size of this region and consequently 
the size of these negative peaks. 



This subsection concludes the study of the perturbative 
caustics with the analysis of the position of cusps in these 
caustics. The position of cusps can be important in sev- 
eral studies such as microlensing itself. In fact cusps are 
surrounded by a region with an amplification even higher 
than that of fold singularities. They also define the exten- 
sion and the shape of the caustic. 

Cusps are defined as the points where the tangent vec- 
tor of the caustic vanishes. In order to find them we must 
set 

2/i {^) = 
2/2 (^) = 

and resolve this system of equations for ?9. 

Let's start with the central caustic. Eqs. ([l^ ) after sev- 
eral steps become: 



(18) 
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" ( A°2 COS 1?- A"i sin i?) (A^^ cos iJ+A^j sin i?) 

(19) 



8 



Valerio Bozza: Perturbative analysis in planetary gravitational lensing 



Jl 



0.0? 



.0 



, 1 : 



-1. 4 -l.i -1 -0. S -0.6 -0.4 -O.t 



-0.0 



0.0? 



0.1 



we have 
mi 



mi 9r 
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Multiplying by r, we get an equation in 
easier to handle: 

— r'' (cos^ d - sm^ i?) H ^^r-- smi?cosi? - 



Inserting (^l]) and solving, we have: 

2 2 
on the higher branch, and: 



± - arcsin 
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2mi 



•d — n ± — arcsin 
2 



Fig. 6. For planets internal to the star's Einstein ring, the 
planetary caustics move on this curve. 

These can simultaneously vanish only if 
de ^ (AOj cos I? - A^i sin i?) ( A^i cos i9 + A% sin d) 
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— ± - arcsm 
2 2 
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(23) 
which is 

+ m2 = 
(24) 

(25) 

(26a) 
(26b) 
(26c) 
(26d) 



Explicitly, this equation is: 



E 



(20) 



sin ?? cos 



on the lower. 

= For external planets, the higher branch is complete 
while the lower is absent, so only the four cusps on the 
higher branch are actually present. For internal planets. 



the two branches are real only near d 



and ■§ 



So the higher branch has the two cusps dX ■& ~ 



3ZL 

2 ■ 

and 



d — 4f , while in the lower one the four cusps (26c) and 



+ 3A0iA«, [(A^i)^ + (AO^)'] (cos^ ^ ~ sin^ ^) + 
A^iV^ [(AOi)' - 3 (4sin3 d - 3 sin ^) + 



+ A^2V^ 



(A02)'-3(A0i)'] (4cos3^-3cos^)} 



(21) 



which, in despite of its cumbersome aspect, can be exactly 
solved in the case of the single planet where it yields the 
four solutions: 



•d — ± arccos 



3mi + 3x1 - ^9ml - Umixl + 9xj, 
Ay'miXp 



(22) 



For planetary caustics, we can proceed in a similar way. 
Multiplying (|l^a) by sini9, (^8|b) by cosi? and subtracting. 



( 26d ) are real and the others are imaginary. Summing up 
we have six cusps distributed in such a way as to form two 
triangular caustics. 



4. Microlensing 

Among the numerous forms of gravitational lensing, mi- 
crolensing is surely one of the most relevant since it opens 
the possibility of probing the galactic structure through a 
directly gravitational investigation. 

Microlensing occurs when the images of a given source, 
produced by a small lens, are too close (typically less than 
10"'^ arcsecs) to be separated by our telescopes. As we can- 
not see but a point image of the source, the only way to 
notice a lensing effect is through a variation of the total 
light flux coming from the observed source. For a point 
lens mass, this variation was found by Paczynski (1986) 
who first thought of galactic microlensing as a new ob- 
servable astronomical phenomenon. For a planetary sys- 
tem the anomalies in amplification patterns do not enjoy 
a full analytical description. Our aim is to use a pertur- 
bative approach to solve this problem and find analytical 
light curves for stars accompanied by their planets. 
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Fig. 7. Typical microlensing curve for a point lens mass. 
J^.l. Paczynski's curve 

Before considering the problem of planetary microlensing, 
it is useful to review the steps to be followed in order 
to get amplification light curves in the event of a single 
mass (Schneider, Ehlers & Falco, 1992). This will help us 
in fixing the problems to be faced. In this case, the lens 
equation takes a very simple aspect: 



y = X 



nil 



(27) 



The total amplification is found by summing the am- 
plification of all images. So the first step is to find these 
images, i.e. the lens equation is to be inverted. Here the 
task is quite easy, because (|^ reduces to a second degree 
equation whose solutions are: 



y 

2 

y 

2 



^ 4mi + 


\y\' 


|y| 


^J 4mi + 


\y\' 



(28a) 



(28b) 



The next step is to compute the amplifications corre- 
sponding to each of these images. According to (j|), these 
are: 



1 



|det J(I±)| 



|i± 



(29) 



It is interesting to study the properties of the two 
images to discover their physical essence (Blandford & 
Narayan, 1986). The image 1+ has positive parity; in the 

limit of vanishing lensing effect ^|y|^ S> 1+ tends to 

y and its amplification becomes unitary. Thus reduces 
to the usual image of the source in the absence of lensing. 
In what follows I'll refer to it as the principal image. I~ 
has negative parity and in the limit of low lensing goes as 
'^i/ |y|: while its amplification is always fj,j- = iJ,j+ — 1. 



I shall call it secondary image as it disappears when the 
lensing effect is not present. Both images are aligned on 
the line connecting the source and the lens: the principal 
image is always external to the Einstein ring, while the 
secondary one is internal to it. 

Now we have to sum up the two amplifications to ob- 
tain the so - called amplification map: 



2rni 



(30) 



4m 1 



This function tells us the amplification corresponding 
to any given position of the source relatively to the lensing 
object. Of course it only depends on the distance because 
of the symmetry of the lens. 

The final step is to make the source move along a rec- 
tilinear trajectory to obtain the complete light curves cor- 
responding to the passage of a massive lens near the line of 
sight of the source (obviously it makes no difference who 
is moving, what counts is only the relative motion). The 
distance |y| is: 



62 



'i2 



(31) 



where b is the impact parameter (the closest approach 
distance) and v± is the projection of the relative speed in 
a plane orthogonal to the line of sight. 

A typical light curve is shown in fig. 0. The height 
of the maximum is found by substituting the impact pa- 
rameter in (|30|). It becomes infinite as b tends to zero. 
Real sources have finite extensions implying integration 
processes smoothing the peak of the curve (Witt & Mao, 
1994). This cut - off becomes evident when |y| is compa- 
rable to the source radius. 

4-.2. The problem of planetary microlensing 

In principle, the procedure for attaining microlensing light 
curves for multiple lenses is the same just expounded for 
a single point lens. First we should invert the lens appli- 
cation, second we have to compute the amplification of all 
the images, then sum up to have the amplification map 
and finally introduce the motion of the source relatively 
to the lens system. But if wc write the lens equation for a 
star with just one planet placed in Xp = {xp, 0): 



mi 



1712 



(X - Xp) 



(32) 



we at once see that the inversion is not possible since one 
must surrender at a fifth degree equation which doesn't 
allow to find the images produced by such a lens. 

A glance at the numerical results can indicate us which 
way is to be taken in the inversion of the lens application 
(p2[). When the source is outside the caustics, only three 
images are present (see fig. One of them is outside all 
critical curves and approaches the source when the latter 
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Fig. 8. In presence of a planet, a source (here represented 
by a four-cornered star) placed outside the caustics gen- 
erates three images (the three little circles in the figure). 



is far enough from the lens. This is indeed the principal 
image. Another image is inside the star's critical curve. It 
is easy to understand that when the mass of the planet 
vanishes this image becomes the star's secondary image. 
The last image is near the planet (inside the planetary 
critical curve when the planet is external to the star's 
Einstein ring). I shall refer to this as the planetary image. 
It is clear that the presence of the planet slightly perturbs 
the principal and secondary image of the star, so that 
their position can be found applying perturbation theory 
to Paczynski's images. The planetary image is completely 
perturbative, since it is not present in the zcroth order sit- 
uation in which the planet is absent, and must be treated 
separately. When the source threads a caustic, two new 
images are formed with opposite parities whose effects are 
similar to those of the planetary image. 

So, the perturbative analysis is likely to be the key to 
solve the problem of planetary microlensing. In the fol- 
lowing two subsections I will use it to discover the images 
and their amplification. Finally, I will build amplification 
light curves and compare them with their numerical coun- 
terparts. 



4-3. Principal and secondary image 



Paczynski's images (^8[) are the starting point for our ex- 
pansion and will be generically indicated by the symbol 
I^"). Let's write the position of the images to the first or- 
der in TO2 as the sum of Paczynski's image and a small 
perturbation AI: 



Id) =l(o)+AI 



(33) 



With this position, in the lens equation expanded to the 
first order in m2 



- r(i) 



mill 



(1) 



r(0) 



r{0) 



(34a) 



y2 



mil. 



(1) 



m2l. 



(0) 



(^! 



(0) 



(34b) 



the planetary term no longer contains the perturbation 
AI. Putting: 



7712 



AO) 



AO) 



m2l. 



(0) 



1 J^p 



(35a) 
(35b) 



and bringing these terms to the left members, we re- 
gain the structure of the Schwarzschild lens equation (27) 
in the variable 1^) for the source position y + Ay. The 
planetary induced perturbation can be thus read as a shift 
in the source position. 1^) has the same expression as I^°) 
evaluated in y + Ay instead of y. The perturbation AI is 
found by expanding 1^) to the first order in Ay: 



A/? 



At/i 



1 ± 



^47771 




2 V 

{yi^yi + 2/2Aj/2) 

^^mi^yi + vl^{yl + ylf 



(36a) 



A/± = 



A2/2 



1± 



V47771 +yi + : 



vl 



2 Vvl 

2?77l7/2 (2/1 Ayi + 7/2 Ay2; 



v/4777i+7/2+y2^(y2^y2)3 



(36b) 
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The upper signs stand for the principal image and the 
lower for the secondary. The expansion parameter m2 ap- 
pears through Ay. 

Now the position of the principal and secondary image 
are known. The most delicate operation is done and the 
door to the planetary microlensing is open at last. What 
remains is only mechanic computation without any con- 
ceptual difficulties. 

The amplification of each image is found by the general 
formula expanded to the first order in m2 (I drop the 
zero from l'"' to simplify notation): 



1 - 



(/1A/1+/2A/2) 



2miTO2 



1 - 



[II - II) ((/i - Xj,f - /|) + 4/1 (Ji - Xp) I\ 

{ii+iif{{h-x,f+iif 



(37) 



This is the sought formula for the amplification of 
the images. Paczynski's amplification multiplies the main 
brackets containing the sum of all perturbations following 
the zero order solution represented by the unity. Two kinds 
of perturbations can be recognized: the first is caused by 
the previously found shift in the image positions AI; the 
second is the consequence of the change of the function 
det J produced by the presence of the planetary term in 
the lens equation. I have dropped the modulus from the 
main brackets because its content is always positive since 
the perturbations are smaller than unity (except for the 
zones where perturbative method is no longer valid). 

As usual, the validity of perturbation theory is limited 
to the regions where perturbations are enough small to 
make sense. So it is necessary a careful examination of the 
denominators of all perturbative terms. The shift terms 
present the distance of the zeroth order image from the 
origin raised to the sixth power. There's no problem for 
the principal image which is always far beyond the Ein- 
stein ring, but this is not true for the secondary image. 
However the "failure" rises in the limit of vanishing lens- 
ing where the amplification of this image is so low to be 
totally masked by the amplification of the principal image. 
When the amplification of the secondary image begins to 
become important, the distance from the origin is largely 
sufficient to eliminate all the problems and have fine per- 
turbations. The shift AI becomes infinite when the source 
passes through the origin; so the region very near the ori- 
gin is the first to avoid. The displacement Ay diverges 
when the zeroth order image approaches the planet as 
could easily be foreseen for a first order perturbation the- 
ory. As regards the terms coming from the alteration of 



det J, there's nothing new; the prescriptions are the same 
as those from the other terms. 

In sum we are allowed to use these amplification for- 
mulae for all source positions being not too near the origin 
or generating images too close to the planet. This hardly 
happens when the source is internal to the caustics. We'll 
see that very reliable results can be obtained up to very 
short distances from the centres of the caustics. 

J^.^. Planetary image 

As previously announced, in this subsection I shall deal 
with the third image. The presence of this image is ab- 
solutely tied to that of the planet. Anyway, Paczynski's 
images can still provide a good starting point for our anal- 
ysis. In fact, if the planet is very far from the star, it too 
will behave as a single lens. In this case, the planetary 
image is nothing but the secondary Paczynski's image for 
a very low mass. In this limit, its distance from the planet 
is of order 7712- So, in our perturbative expansion, we are 
encouraged to search for images with distances from the 
planet of order m2 . Let V be the position of the planetary 
image. We have: 



AF 



(38) 



with AP of order 7712. Saving only the lowest order, the 
lens application reads: 



mi 

2/1 = Xp 



maA/f 



2/2 



X, {Mlf + (A/f )^ 
maA/f 



(A/f ) + (A/f 



(39a) 
(39b) 



These equations can easily be solved. The solution is: 
_ TO2 (2/1 - 2/p) 



(2/1 - Vp) + 2/1 

"722/2 



(2/1 - 2/p) + vl 



(40) 
(41) 



where 2/p = Xp — ^ is the zeroth order position of the 
planetary caustic already rising in former discussions. 

As ever the amplification is calculated by expanding 
(0). The lowest order result is: 



(2/1 - 2/p) + 2/1 



(42) 



Notice how this formula is much more simple than 
other images amplification. 

The denominators in these expressions vanish when 
y — > yp. Consequently the perturbative method fails when 
the source is very close to the centre of the planetary caus- 
tic. 
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4-. 5. Perturbative light curves in planetary microlensing 

Once we have found the amplification for each image, in 
order to obtain the microlensing amplification map we 
must sum up the components coming from the three im- 
ages. However, we see that the contribution to the total 
amplification of the planetary image is of the second order 
in m2. Since we are only considering first order corrections 
to Paczynski's curve, this contribution is to be ignored. 
Therefore, from now on, we shall confine ourselves to the 
principal and secondary images only. 

One consideration is for the two hidden images coming 
out when the source crosses a caustic. If the event regards 
the planetary caustic, the two images can be found by 
carrying further the expansion ( |38| ) . The new images arise 
from higher order solutions and their amplifications will 
also be of higher orders. So we don't worry about them. On 
the contrary, if the source crosses the central caustic, the 
new images appear near the star's Einstein ring, far from 
any possible starting point for a perturbative expansion. 
As we are not taking them into account, we cannot expect 
to obtain good results inside the central caustic. Anyway, 
central caustic crossing events are very improbable, since 
the extension of this caustic is ml/mf times the star's 
Einstein ring. 

Building light curves presents no difficulty. Chosen one 
source trajectory, it suffices to parameterize yi and t/2 in 
the amplification map properly. This is no longer a func- 
tion of the radial coordinate because there is no more ro- 
tational symmetry. 

To account for the finite size of the source a simple nu- 
merical integration of the perturbative amplification map 
on the source area at each point of the trajectory can be 
performed. The curves thus obtained can be compared to 
numerical ones given by "inverse ray shooting" algorithm. 

All the results I show in this paper regard a system 
constituted by a star with mass mi = 1 and a Jovian 
planet (m2 — 10^'^). This choice has been made in or- 
der to put in better evidence planetary perturbations and 
to test the perturbative approach in the least favourable 
situation. Obviously with Earth - like planets things can 
only go better. 

Let's start with an external planet. In fig. ^ the planet 
is in Xp = 1.2. The trajectory chosen for this first test is 
shown in fig. |9^ and has impact parameter 0.5. The nu- 
merically attained light curve is displayed in fig. ^Jd. The 
source used for this curve has radius 0.045. In a standard 
observation towards the bulge of the galaxy {Dq l ^ 8kpc, 
Dos lOkpc), this value would correspond to a giant 
about 43 times greater than the sun. Here the presence 
of the planet is responsible for the little peak on the left 
of the maximum of the curve. Fig. ^ represents the per- 
turbative light curve for a point source moving along the 
same trajectory. If we perform the numerical integration 
of the perturbative amplification map, as previously said. 



the perturbative light curve |9|d becomes indistinguishable 
from the numerical one. 

This is a very encouraging result, so let's choose other 
trajectories to see other tests. In fig. |l^the position of the 
planet is the same but the trajectory passes between the 
planetary caustic and the central caustic at a minimum 
distance of 0.2. The peak in the numerical curve is 
very close to the maximum. The point source perturbative 
curve |lO|c presents a sharp peak which assumes the right 
proportions after the integration in [lO|d. 

At this point, let's see what happens when the source 
crosses the planetary caustic. In fig. pl]a, the impact pa- 
rameter 0.4 allows the crossing. The peak in the numerical 
curve |l^b becomes considerably high. In the point source 
perturbative curve pT| c the peak is very sharp (it would 
diverge at the centre of the caustic yp). However, the in- 
tegration over the source surface still succeeds in reporting 
this peak to the right size and shape. 

Now, let's consider an internal planet {xp — 0.8). The 
region between the couple of planetary caustic is charac- 
teristic for its high de-amplification. This produces nega- 
tive peaks on light curves such as the one shown in fig. 
p2[ b corresponding to the trajectory in fig. p^a. It is in- 
teresting to see that the perturbative method reproduces 
this situation with the same great accuracy proved in the 
former situations. As ever, the point source peak in [l^ is 
smoothed by finite source effect in |l2|d. 

In fig. ^ the impact parameter is 0.25 and things go 
perfectly as previously. 

Finally, let's consider caustic crossing in this case. Fig. 
shows a trajectory very close to the planetary caus- 
tics. The "inverse ray shooting" curve presents a large 
de-amplification preceded and followed by little positive 
peaks. The perturbative curve |ljc is characterized by the 
same situation but the de-amplification is so high to make 
the total amplification become (unphysically) negative. 
Now let's see what happens with a finite source. Because of 
its extension, part of the source hits the centre of the caus- 
tic Yp where the perturbative amplification map diverges. 
This is a hard problem for the numerical integration which 
becomes very unstable in this zone, so the bottom of the 
de-amplification region of the light curve |ljd cannot be 
taken as significant. However, we see that things go fairly 
well even in this extreme situation. 

5. Conclusions 

The success of perturbative theory in planetary lensing 
cannot but impress for the simplicity of the calculations 
involved and the surprising accuracy of the results even in 
the hardest situations. 

In the derivation of the caustics of a planetary sys- 
tem, by a simple idea and very few passages the complete 
structure of these curves has been easily obtained. The al- 
most complete insensitivity of the perturbative approach 
to the number of the planets allows complete descriptions 
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Fig. 9. Besides a star with unitary 
mass placed in the origin, here is a 
Jovian planet (m2 = 10 ~'^) in Xp = 
(1.2; 0). The trajectory of the source 
shown in (a) has impact parameter 
0.5. (b) is a numerical light curve 
obtained by "inverse ray shooting" 
for a source 43 times greater than 
the sun. (c) is the perturbative light 
curve for a point source and (d) is 
the same curve after a numerical in- 
tegration over the source extension. 
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Fig. 11. The lens system is the same 
again but the impact parameter 0.4 
makes the source cross the planetary 
caaustic. 
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Fig. 12. Now the planet is in Xp = 
(0.8,0). The three light curves are 
obtained in the same ways described 
in fig. H The impact parameter is 
0.6. 
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2 Fig. 13. Same lens system as in the 
previous figure. The impact param- 
eter is 0.25. 
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Fig. 14. Caustic crossing with im- 
pact parameter 0.4. 
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of planetary systems without any loss of generality. Also 
many important physical assertions can be stated thanks 
to these results. The fact that the shape of the central 
caustic is largely given by a linear superposition of the 
effects of the single planets is indeed remarkable. 

In planetary microlensing the results are even exalting. 
The perturbative amplification map allows the construc- 
tion of very fine light curves. In the derivation of the am- 
plification map I have dealt with only one planet for the 
sake of simplicity. Yet the generalization to an arbitrary 
number of planets is immediate because in the first order 
domain a superposition principle is here valid as well. For 
point sources, light curves can be attained in a completely 
analytical way, while for finite sources I have resorted to 
numerical integration until now. Work to englobe finite 
source effect in the analytical description is in progress. 
When these curves are available, the extraction of param- 
eters of planetary systems from microlensing light curves 
will start on more solid analytical bases. Also it could be 
possible to use the analytical expressions in experimental 
fits, though the large number of parameters would greatly 
affect the uncertainties in their determination. 
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